Panorama of Nodal Superconductors 

Hyekyung Won/ Hyunja Jang/ David Parker/ Stephan Haas/ and Kazumi Maki^ 

^Department of Physics, Hallym University, Chuncheon 200-702, South Korea 
^Department of Physics and Astronomy, 
University of Southern California, Los Angeles, CA 90089-0484 USA 

(Dated: February 2, 2008) 

Abstract 

Since 1979, many new classes of superconductors have been discovered, including heavy-fermion 
compounds, organic conductors, high-Tc cuprates, and Sr2Ru04. Most of these superconductors 
are unconventional and/or nodal. Therefore it is of central importance to determine the symmetry 
of the order parameter in each of these superconductors. In particular, the angular-controlled 
thermal conductivity in the vortex state provides a unique means of investigating the nodal struc- 
ture of the superconducting energy gap when high-quality single crystals in the extremely clean 
limit are available. Using this method, Izawa et al have recently succeeded in identifying the en- 
ergy gap symmetry of superconductivity in Sr2Ru04, CeCoIns, k-(ET)2Cu(NCS)2, YNi2B2C, and 
PrOs4Sbi2. 

PACS numbers: 



1. Introduction 

Since tlie appearance of anisotropic superconductors the determination of their gap sym- 
metries has been one of the central issues. For the high-Tc cuprate superconductors 
angular-resolved photoemission spectroscopy (ARPES)P 



and Josephson interferometry 



provide a definitive signature of dj,.2_y2-wave superconductivity. However, these meth- 
ods have not yet been extended beyond the high-Tc cuprates. In a remarkable paper 
Volovik[3] showed that the quasiparticle density of states in nodal superconductors in a 
magnetic field is calculable within a semi-classical approximation. Here the Doppler shift ^ 
due to the supercurrent circulating around vortices plays a crucial role. The resulting \/H- 
dependence of the s pec ific heat in the vortex state has been observed experimentally in 
YBCoQ,Q, LSColnj, fi:-(ET)2 salts 12|, and Sr2Ru04Q |l4|. Here H is the magnetic 
field strength. 

TM. app.oacK ha. been expended i. .an. d.ee.io.. K.be. and Hi.chMd0 eonsid- 

ered the effect of non-zero temperature and established the scaling law, as formulated by 
Simon and LeeQ]. In particular, Kiibert and Hirschfeld found the scaling function within 
the semi- extended this approach for the superfiuid 

density, the spin susceptibility, and the nuclear spin lattice relaxation rate. Measurements 
of the thermal conductivity tensor in the vortex state provide an advantage, since the heat 
current direction provides further directional information. Kiibert and Hirschfeld 3] and 



Vehkter et al 



17| have performed analyses in these directions. However, these authors did 



not take the spatial average over the vortex lattice at the same time as the average over 
the Fermi surface, but instead assumed a local thermal conductivity ^^(r). It is known 
that the local quasiparticle density of states mostly consists of bound states around vortex 
cores. These bound states do not contribute to thermal conductivity when H/Hc2 <^ 1. 
On the other hand, the nodal excitations which do contribute to the thermal conductivity 
run through many unit cells of the vortex lattice. Furthermore Vehkter et al chose a rather 
specific circular Fermi surface instead of the cylindrical Fermi surface commonly used in 
modeling the high-Tc cuprates. These problems were addressed in references 18-21. In these 
works simple expressions for the thermal conductivity are found in the superclean limit 
((FA)^/^ <^ v\J eH <^ A(0), where F is the quasiparticle scattering rate in the normal state 
and V is the effective Fermi velocity). 

Intuitively we understand that the Doppler shift generates quasiparticles in a plane per- 
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pendicular to the vortex axis and the field direction. When this plane intersects the nodal 
directions, there is enhancement of the available quasiparticles, which we call nodal excita- 
tions. Therefore, by changing the field direction one can sweep the Fermi surface with the 
plane associated with the Doppler shift. As this plane meets the nodal directions both the 

nnn 

specific heat and thermal conductivity are enhanced [18|, [191, |20||. 

In particular when the Fermi surface is quasi-2D or cylindrical, rotating the magnetic 
field within the conducting plane yields valuable information regarding the gap symmetry. 
This was shown clearly in early experiments on YBCoQ, Q, 25 1 and more recent experi- 
ments on Sr2Ru04j3|, CeCoInsj^^l, and /t-(ET)oCu(NCS)o j28|. Such experiments, probing 
the angular dependence of the superconducting order parameter, indicate f-wave supercon- 
ductivity for Sr2Ru04, d3,2„j^2-wave superconductivity for CeCoIus and k-ET2Cu(NCS)2, 
and s+g-wave superconductivity for YNi2B2C and PrOs4Sbi2. These order parameters are 
shown in Fig. 1. For A(k) of PrOs4Sbi2 see Fig. 10. 






FIG. 1: 2D f-wave, dx2_y2-wave, and s+g-wave symmetry. 

In the following sections we focus on the salient features of superconductivity in high-T^ 
cuprates, Sr2Ru04, k— (ET)2 salts, YNi2B2C, and PrOs4Sbi2. 
2. High-Tc Cuprate Superconductivity in a Nutshell 

The discovery of high-Tc cuprate superconductivity in LBCO by Bednorz and Muller[ 
in 1986 took the superconductivity community by surprise. The subsequent excitement and 



confusion are well documented in a textbook by Enzj3o|. In 1987 P.W. Anderson 31 1 pro- 
posed the 2-dimensional one-band Hubbard model and his famous dogma. His central theme 



is to understand the superconductivity in the presence of the strong Coulomb repulsion. In 
the meantime (ia;2_y2-wave symmetry of botli t_he Jiole and the electron-doped high-Tc cuprate 



superconductivity has been established ^, 



. Furthermore, the mean- fie 
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d theory (i.e. 



A 
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so within 
have 



the generalized BCS theory) for d-wave superconductivity works well 
the framework of the BCS theory of d-wave superconductivity, May Chiao et al 
derived the crucial parameter A{0) /Ep of optimally doped YBCO and Bi-2212 through the 
measurement of thermal conductivity at T < IK. Here A(0) is the maximum value of the 
energy gap at T = OK and Ep is the Fermi energy. They found A{0)/Ep ~ 1/14 and 1/10 
for YBCO and Bi-2212 respectively. Regrettably there are no experimental data available 
indicating the ratio A{0) /Ep in the underdoped and overdoped region of YBCO and Bi- 
2212, though we have no reason to worry that this ratio becomes substantially different. 
First of all, these values tell us that high-Tc cuprate superconductivity is very far away 
from the Bose-Einstein (BE) condensate, but is within the BCS regime. The BE condensate 
clearly requires A(0) ~ Ep. Secondly, making use of Ginzburg's criterion the fluctuation 
effect in high-Tc cuprate superconductivity should be at most of the order of a few percent. 

Thirdly, these A{0) /Ep values are incompatible with the assumption that A(0) ~ Ep, 
which was made in solving the Bogoliubov-de Gennes equation in the vortex state of d-wave 



superconductors 
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38, 



39|. In particular the approximation A(0) ~ Ep appears to knock 



off all the bound states around a single vortex in d-wave superconductors. On the other 
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4l|. 



hand, for A{0) /Ep = ^, one finds hundreds of bound states around a single vortex 

In fact, the local density of states around a single vortex looks very similar to the one 
found for s-wave superconductors j^]. Second the approximation A(0) ~ Ep appears to 
increase the y/H term in the specific heat by a factor of 10 to 30 j^]. Therefore we stress 
that if we limit ourselves to the region T ^ (|v ■ q|) ^ ^(0), the quasiclassical approach 
as discussed in Refs. 18-20 provides the most reliable result so far available. 

Unfortunately the approximation A(0) ~ Ep is very popular among the superconductiv- 
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ity community because it makes the calculation much simpler 

As to the observation of these bound states around vortices in YBCO and Bi2212, only 
a few bound states, if any, are observed^, |42|- But this appears to be due to the fact 



or 



that the vortex core is filled with other order parameters like antiferromagnetismj4 
charge density wave. 

In a popular paper Bob Laughlinjs^ proposed a very intuitive picture of superconduc 
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tivity in the presence of the Mott insulator. We think that another interesti ng q uestion is 
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53|. 



how the superconductivity can survive in the presence of d-wave density wave|51l. 

In section 4 we shall discuss a very similar conflict between two order parameters in the 
organic superconductor k-(BEDT-TTF)2X. 
3. Superconductivity in Sr2Ru04 

Superconductivity in Sr2Ru04 was discovered in 1994[5^. The surprising prediction of 



triplet p-wave pairingand related chiral symmetry breaking 



rotation experiments 



55| was verified 



jy muon spm 



56l | and a flat Knight shift as seen by NMR measurements 



triplet superconductivity implies clapping collective modes 



571]. Also, the 
and half-quantum vortices [5^ 



as topological defects. Since sample quality has improved and single crystal Sr2Ru04 with 
Tc — 1.5 K have become available, experiments have found clearly nodal structures jisl. IgoIi . 
This is inconsistent with the initiall y p roposed fully gapped p-wave model. Consequently, 



several f-wave models were proposed 61, 62, 63 1 . 

As shown in Ref. 19, most of these models are consistent with the specific heat datallSf 

. The p-wave model cannot account for these 
measurements. Furthermore, ultrasonic attenuation data by Lupien et al|sJ] eliminates 
models with vertical nodal lines parallel to k^. This leaves only the 2D f-wave model with 
energy gap 



A(k) = Ae±*^cos(cA; 



(1) 



where e^* 

Zhitomirsky and Rice|65l. 



have proposed a multi-gap model where one of the supercon- 



ducting order parameters associated to the a and (3 band has horizontal nodes at the Bril- 
louin points j67|. However, this multi-gap model gives a two-fold-symmetric angular depen- 
dence (~ cos(20), where (p is the angle between the heat current and the magnetic field) ten 
times larger than observed experimentally 26|. Moreover, the two-gap model cannot give uni- 
versal heat conduction as observed in k ^^ | 6q . More recently, further tests of the 2D f-wave 



model for Sr2Ru04 were proposed[69|,l70|,l71|. If the f-wave model is established for Sr2Ru04, 
it implies that the 2D model used in the high-Tc cuprates is inapplicable to Sr^RuQ4. From 
the point of view of the electronic interaction Sr2Ru04 should be a 3D system[22]- Also it is 
possible that p-wave pairing is forbidden in the electronic systems. It appears we now have 
3 f-wave superconductors: UPtsI?^, Sr2Ru04, and UNi2Al3[2^. Finally, it is important to 



identify the gap symmetry of the organic superconductor (TMTSF)2PF6 for which triplet 
pairing has aheady been estabhshed[75|. 

Very recently Deguchi et al 3] have reported specific heat data for 0.12 K < T < 0.51 
K and for magnetic field 0.05 T < B < 1.7 T and rotating within the a-b plane. One of the 
most surprising data is the fourfold term with cusps at = 0, 7r/2, etc., which appears only 
for B = 0.30r,0.60r and O.OOT at T = 0.12K. For data at T = 0.31K the cusp feature 
becomes less clear. 

Perhaps the most interesting question is 1) Does this include 2 energy gaps?, and if so, 
2) What is the gap symmetry of this new gap? 

If we compare their Fig. 3 for T = 0.12K, the angular dependence is very similar to 
what one expects in s+g-wave superconductivity in the presence of impurities |77t|. In s+g- 
wave superconductors the impurity scattering induces a small energy gap. Therefore, the 
most natural interpretation of the above data is that the specific heat has picked up four 
point-like mini-gaps located at 6 = 7i/2 and = 0,7r/2, etc. The size of the minigap 
would be around 0.3 K. Then it is rather difficult to accommodate this with the Miyake- 
Narikiyo model|3]. Also this quasi- nodal structure is difficult to accommodate with the 
f-wave order parameter above. Perhaps this experiment indicates the presence of the second 
energy gap? The gap function A(k) with minigap is shown in Fig. 2. Here we take 
|A(k)| ~ |cos(x)|(l - 1.8cos(40) cos2(x) + 0.81 cos^(x))^/^. 




FIG. 2: Proposed Gap Function for Sr2Ru04 
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4. Gossamer Superconductivity in k-(BEDT-TTF)2X? 

The organic superconductors fi:-(BEDT-TTF)2X with X = Cu[N(CN)2]Br, Cu[N(CN)2]Cl 
and Cu(NCS)2 have the high est superconducting transition temperature = 10-13 K 
among organic conductors 80| . There are many parallels between high-T^ cuprates and k- 
(BEDT-TTF)2X; the quasi-two dimensionality of the Fermi surface, and the proximity to the 



antiferromagnetic state. More recently angular dependent STM 



8l| and angular dependent 



thermal conductivity measurements in the vortex state j28| indicate dx'^-y'^ superconductiv- 



ity. The nodal lines runl28ll in the diagonal direction of the b-c (i.e the conducting) plane. 

nn 

Although d-wave superconductivity has been speculated theoretically the diagonal 
lines come as a surprise. This indicates that perhaps the exchange of an antiparamagnon is 
not adequate to generate d-wave super conduct ivitvjs^ . 

Indeed from the thermal conductivity data from fi;-(BEDT-TTF)2Cu(NCS)4, we can de- 
duce that A(k) = A(cos(20) —0.067). This d-|-s-wave symmetry is somewhat similar to that 



of YBCO 



85|. 



Perhaps more surprising is the sensiti vity of the superconductivity to the cooling rate, 
which has been studied by Pinteric et aljS^ in K-(BEDT-TTF)2Cu[N(CN)2]Br. To make 
the story simple we consider two extreme cases only. We define the relaxed sample as the 
one cooled very slowly down to 10 K. For example this sample is left for 3 days in liquid 
N2 and then cooled down for a fraction of a degree K per hour down to 10 K. The other 
extreme is the quenched sample, wherein the sample is dropped in liquid N2, then cooled 
down to 10 K within a few hours. 

Of interest, the superconducting transition temperature Tc shows little dependence on the 
cooling procedure, with Tc varying by only a few percent. More surprising is the sensitivity 
to the cooling rate of the superfiuid density as measured by magnetic penetration depth. 
The superfiuid density of the quenched sample is only 1-2 percent of the relaxed sample. 
Also the temperature dependence of the relaxed sample exhibits a T-linear dependence 
typical of d-wave sup er conductors S|- other hand, the superfiuid density of the 

quenched sample is very flat for T < 0.2 Tc, which may be interpreted as a sign of s-wave 
superconductivity. Until recently there were debates on the symmetry of superconductivity 
in K-(BEDT-TTF)2 salts: d-wave versus s-wave. We therefore wonder if this controversy 
originates from the difference in the cooling rate. For example, Elsinger et aljs^ did not 
describe how their samples were cooled down. 
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Recently the effect of the coohng rate on the normal state of the three k-(BEDT-TTF)2 
salts has been reported 9^] ■ These organic conductors go through a glassy phase when they 
are cooled down from room temperature to liquid N2 temperature (70 - 90 K). 

The BEDT-TTF molecules have ethylene groups attached to each end. For T > lOOK 
these ethylenes are rotating freely. As the temperature decreases below 70 K these ethylenes 
cannot move freely but are settled in their equilibrium configurations. Therefore it is very 
likely that in the relaxed samples these ethylene groups are relatively well ordered while in 
the quenched samples they are oriented randomly. 

But it is not known at present how the randomness of the ethylene groups interferes with 
the quality of the superconductivity. From the insensitivity of Tc to the cooling procedure 
we infer that this cannot be the simple effect of disorder. Also Pinteric et al have shown 
the superconductivity to be homogeneous. More likely is that the disorder controls the 
appearance of another order parameter, a "hidden order parameter", like unconventional 
density wave (UDW). This forces the superconductivity to coexist with UDW. 

As we have discussed in section 2 we call this type of superconductivity "gossamer 
superconductivity" Q^. Unfortunately at present we do not know what kind of order pa- 
rameter is appropriate to characterize the glassy phase, although we suspect it should be 
UCDW or USDW. We believe that this is the most fascinating question in k- (BEDT-TTF) 2 
salts. 

5. Superconductivity in YNi2B2C 

Superconductivity in YNi2B2C and LuNi2B2C was discovered in 1994 Recent interest 
has focused on these materials because of their relatively high superconducting transition 
temperatures of 15.5 K and 16.5 K, respectively. There is a substantial s-wave component 
in the order parameter of these compounds. This was shown by substituting Ni with a 
small amount of Pt and observing the subsequent opening of the quasiparticle energy gap as 
observed in specific heat measurements |93l| . Nevertheless, recent experiments clearly indicate 
the presence of nodal excitations 

0,0- 

Furthermore, the upper critical field within the a-b plane exhibits a clear four-fold 



symmetry 
pears to be 
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The simplest model to describe the superconductivity in YNi2B2C ap- 
A(k) = ^A(l - sin^ e cos(40)) (2) 
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or s+g-wave superconductivity. Here 6 and are the polar and azimuthal angles, respec- 
tively, describing the direction of k. The corresponding A(k) is shown in Fig. 1(c). The 
quasiparticle density of states is given by 



G{E) = N{E)/No = \E\Re{ 



(3) 



where ( ) means (l/47r) J dVt 

\E\/A < 1, we obtain 



(^2- A2(k))V 

. This is shown in Fig. 3. In particular, for 



^ ' 4A 16 ^ ^ ^ 



(4) 
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FIG. 3: Quasiparticle density of states of an s+g-wave superconductor. 
This gives rise to power laws in the specific heat and other quantities as follows 
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(9) 



9 



where is the Sommerfeld coefficient, and C(3) ~ 1.202 is the Riemann zeta function. We 
note that the presence of point nodes in the a-b plane creates an anisotropic temperature 
dependence in the superfiuid density. 

These specific heats, spin susceptibihty and anisotropic superfiuid densities are compared 
with those of a d-wave superconductor in Figs. 4, 5, and 6(a) and 6(b). 
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FIG. 4: Specific heats normalized by normal state specific heat are shown for s-|-g-wave,s-wave and 
d-wave superconductors. 7iv is the Sommerfeld constant. 
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FIG. 5: The spin susceptibility normalized by that of the normal state are shown for s+g-wave, 
s-wave and d-wave superconductors. 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

t=T/T^ t=T/T^, 

FIG. 6: Superfluid density in the ab-plane, and parallel to the c-axis, for s+g and d-wave super- 
conductors 

In the vortex state and in the superclean limit (F <^ v^/eH, where v = ^/v^ and F is 
the scattering rate in the normal state) the specific heat and other quantities are given by 



where 



and 



XsjT) 
Xn 

PsabjT) 
Psab{0) 



2A 



HO,, 



1 - 



2A 

3vVeH 



4A 



X 



X 



4ln{2/x) 



TT 



A 



J(^,0) = -sin2^sin2 0)V2 + (1 _sin2 0cosV)'^') 



(10) 

(11) 
(12) 
(13) 

(14) 

(15) 



Here 6 and (p are the polar and azimuthal angles describing the direction of the magnetic 
field. 

Recent thermal conductivity [3] and specific heat data 102 1 establish experimentally this 
striking angular dependence. In Fig. 7 we show the experimental data together with the 
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theoretical expression. The cusps at = 0,7r/2, etc. for 9 = tt/2 clearly indicate the 
presence of point nodes. In the presence of line nodes the angular dependence of Kzz is well 
approximated by cos(40). As isseen from the denominator of Eq. 13, the effect of impurity 



scattering is very unusual 



103, 



104|. 




& 45 




T3S 



FIG. 7: Experimental and theoretical angular dependence of I {9, 



Regarding Raman scattering, we present here results from a concurrent work 



lOlj l wherein 



theoretical calculations, based on the s+g symmetry, for th e mo des Aig, Big and B2g are 



performed. The comparison between theory and experiment |105l | is shown in Figure 8, with 
the theoretical results on the left and the experimental results, taken at 6 K, on the right. 
Strong agreement is seen for the whole energy range in the Aig mode, with fair agreement 
for the other two modes. Notably absent from the -B23 data is the secondary cusp feature 
found at = 2 A; we believe that this may be a temperature- related effect and eagerly await 
the results of lower-temperatur e exp eriments. 

From the experimental data jl05 | we can extract A(0)=50.4 K and 64.7 K for YNi2B2C 
and LuNi2B2C respectively. On the other hand, the weak coupling theory gives A(0) = 42.2 
K and 43.3 K respectively, where we need A(0)/Tc = 2.72. Therefore we may conclude that 
YNi2B2C is close to the weak coupling limit whereas LuNi2B2C is in the moderately strong 
coupling limit. 

In the absence of a magnetic field the quasiparticle spectrum in the presence of impurities 
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FIG. 8: Theoretical Raman spectra for the Ai^, B\g and Big modes using the s+g model at T = 
(left side) and experimental Raman spectra for YNi2B2C taken at T = (right side) are shown. 



is determined by |infi | 



io = io ^ Yio( 



(A - A//2)2 - (cD)2 



(16) 



and 



A = A/2 + r(^^M^) 



(17) 



(A - A//2)2 - 



where u and A are the renormalized Matsubara frequencies and order parameter, respec- 
tively, and ( ) means l/(47r) J dQ . 

Then the quasiparticle density of states is given by 



G{E) = N{E)/No = \E\ Re{ 



1 



(18) 



'Co^ - (A - A//2)2 

where G{E) is evaluated at uj = E. The DOS for a few F are shown in Fig. 8. Here we 
have used the so-called Born limit due to the presence of a substantial s-wave component in 
A(k). The unitary limit gives virtually the same result as the Born limit. The most unusual 
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FIG. 9: Density of states for various impurity scattering levels 
feature one notices is the immediate appearance of an ener gy ga p when T ^ 0. This is in 



contrast to the usual nodal superconductors with line nodes |107l |. In particular, the energy 
gap is given by ujg = r(l + 2r/A)~^. This has a number of consequences Q]. First of all, in 
the absence of a magnetic field both the specific heat and the thermal conductivity vanish 
exponentially. 

C,/T,/s:,,/T~(a;,/T)3/2e-/^-« (19) 

There is no universal heat conduction unlike other nodal superconductors 107 . 108^ . Thi s 
remarkable effect of impurity scattering on nodal excitation is clearly shown by Kamata jlOoj ]. 
In Fig. 9 we show the thermal conductivity data of Y(Nii_2^Pt2:)2B2C with x=0.05. As is 
readily seen the fourfold term typical of the pure YNi2B2C has vanished completely. For 
example, the quasiparticle density of states in the presence of both a magnetic field and 



impurities is given by 



221 



G'(H,r) = xcos-'iy)eil-y) 



(20) 



where x is defined in Eq.(14) and y = 
Further we obtain 

X 



((1 - yy - %cos-\y) - y^l^^mi - y) (21) 



2cosh~^(l/x) 



,x 



2 cosh ^{l/x)^x 



-)\cos~\y') 



l-y'^)e{l-y') 



(22) 
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FIG. 10: Thermal conductivity data for Y(Nii_a;Pti:)2B2C with x=0.05. 



where 



1 dVeH 



X 



1 — sin 9 cos 



(23) 



vr A ' ' Ax' 

These expressions indicate clearly that the nodal excitations are eliminated when y > ^J2. 
6. Puzzle of Superconductivity in PrOs4Sbi2 

New heavy- fermion superc ondu ctivity in the skutterudite PrOs4Sbi2 with T^, = 1.8 K 
was discovered quite recently |llO|. The existence of two distin ct superc onducting phases 
and the presence of point nodes in A(k) are of gre at interes t 111 . 112j |. Indeed, recent 



i3 



indicate the presence of two 



angular-dependent magnetothermal conductivity data 
phases, with 4 point nodes in the ab-plane in phase A and 2 point nodes parallel to the b 
axis in phase B (see Fig. 10). 

In order to accommodate the observed nodal structure within the cubic symmetric crystal, 
the following order parameters are proposed: 



AA(k) 



Afl - kt 



A5(k) = A(l-A;,^) 



(24) 
(25) 



for the A and B phases, respectively. 

From an experimental point of view, it must be acknowledged that whether the super- 
conductivity is the spin singlet or the spin triplet is unclear. Here we have assumed the 
spin singlet symmetry. However, a muon spin rotation experiment indicates the presence 
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FIG. 11: Gap functions for A and B phase (right side) of PrOs4Sbi2 



of remanent magnetization which may be in dicat i ve of the spin triplet symmetry 115l | . In 
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deed alternative models have been proposed 
describe the observed angular dependence of the thermal conductivity. 
The models given in Eq.(24) and (25) have the T^-specific heat 

Cs j 27C(3)(T/A)/47r for A phase 
7ivT I 27C(3)(T/A)/87r for B phase 

Of more interest is the anisotropic superfiuid density 
Ps\\{T) ^ 37rln2T 



118|, although these models cannot 



(26) 



114l |. In particular, in the B-phase 



PswiO) 



i6 a' 



and 



1 - 



TT 



tV 2l7rC(3) (T^ ^ 



+ 



(27) 
(28) 

(29) 



Ps±(0) * 16 128 V^, 

where the suffixes || and _L indicate directions parallel and perpendicular to the nodal direc- 
tions. 

Very recently surprising superfiuid density measurements of PrOs4Sbi2 were 



reported |ll8l |. Chia et al measured the superfiuid density parallel to each of the three 
crystal axes and found that the superfiuid densities are isotropic and decrease like at 
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low temperatures. Both the isotropy and the dependence are incompatible with all mod- 
els for PrOs4Sbi2 proposed to date. However, if one assumes that the nodal points in the 
B-phase are always aligned parallel to H when the samples are cooled in a fixed magnetic 
field, our model for the B-phase fits the experimental data precisely. This is plausible, since 
the magnetic field is the only symmetry-breaking parameter in the present model. Also, a 
simple analysis of the effect of the magnetic field on the superconducting order parameter 
tells us that the effect of the magnetic field is minimized when the field is parallel to the 
nodal directions. 

In the presence of a magn etic field the specific heat and the thermal conductivity in the 



superclean limit are given by 



TlXi/A (30) 
xa 1 + 7rxA/2 + 31xy40 



2ln{2/xA) 1 + 31xy64 



/tf^ _ xb 1 + 7rxB/12 + 31a;y40 
'k^ ~ 2ln{2/xB) l + 31x|/64 



(31) 
(32) 



where 



XA = ^^^^((1 - sm\e) cos2(0))^/2 + (1 - sm\e) sm\(f))y/^) (33) 

and 

XB = ^^{l-sin\9)sin\4>)Y/' (34) 

Here we assumed that A(k) is given by equations (24) and (25) for the A-phase and 
B-p hase, respectively. The above expressions are consistent with the thermal conductivity 
data |ll3l |. although experimental data at lower temperatures would help to further clarify 
this point. 

Therefore the model proposed in Refs. 112 and 113 appears to be the most consistent 
with the available data. On the other hand we have assumed the spin singlet pairing, while 



some experiments appear to indicate spin triplet pairing 115[. Further clarification on this 
question is highly desirable. 
7. Outlook 

In the past decade we have witnessed the identification of the gap symmetryof the high-T^ 
cuprates through ARPES^, 0] and phase sensitive Josephson interferometryP,y]. However, 
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these techniques do not appear to be practical to use on heavy-fermion superconductors or 
organic superconductors. 

In recent years measurements of the angular dependent magnetothermal conductivity 
have provided a unique alternative means to explore the gap symmetry of superconductors 
with Tc ~ 1-2 K. In this way the gap symmetries of Sr2Ru04, CeColns, k—{ET)2Cu{NCS)2, 
YNi2B2C and PrOs4Sbi2 have been identified. The next step will be to interpret these gap 
symmetries in terms of available interaction terms in these systems. 

In this journey we have discovered that superconductivity with mixed representations 
play a crucial role in both YNi2B2C and PrOs4Sbi2. This is very surprising, but it appears 
that one must accept this new development. Clearly this will open up a new vista in the 
rich field of nodal superconductors. 
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